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Abstract 

We study the Lyapunov exponents for a moving, charged particle in a two-dimensional 
Lorentz gas with randomly placed, non-overlapping hard disk scatterers placed in a 
thermostatted electric field, E. The low density values of the Lyapunov exponents 
have been calculated with the use of an extended Lorentz-Boltzmann equation. In this 
paper we develop a method to extend these results to higher density, using the BBGKY 
hierarchy equations and extending them to include the additional variables needed for 
calculation of Lyapunov exponents. We then consider the effects of correlated collision 
sequences, due to the so-called ring events, on the Lyapunov exponents. For small 

1 



values of the applied electric field, the ring terms lead to non-analytic, field dependent, 
contributions to both the positive and negative Lyapunov exponents which are of the 
form e 2 lne, where e is a dimensionless parameter proportional to the strength of the 
applied field. We show that these non-analytic terms can be understood as result- 
ing from the change in the collision frequency from its equilibrium value, due to the 
presence of the thermostatted field, and that the collision frequency also contains such 
non-analytic terms. 

KEYWORDS : Lyapunov exponents; Lorentz gas; Extended Lorentz-Boltzmann 
equation; BBGKY hierarchy equations; Long time tail effect. 



1 Introduction 

The Lorentz gas has proved to be a useful model for studying the relations between dynamical 
systems theory and non-equilibrium properties of many body systems. This model consists of 
a set of scatterers that are fixed in space together with moving particles that collide with the 
scatterers. Here we consider the version of the model in two dimensions where the scatterers 
are fixed hard disks, placed at random in the plane without overlapping. Each of the moving 
particles is a point particle with a mass and a charge, and is subjected to an external, uniform 
electric field as well as a Gaussian thermostat which is designed to keep the kinetic energy of 
the moving particle at a constant value. The particles make elastic, specular collisions with 
the scatterers, but do not interact with each other. The interest in the Lorentz gas model 
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stems from the fact that its chaotic properties can be analyzed in some detail, at least if 
the scatterers form a sufficiently dilute, quenched gas, so that the average distance between 
scatterers is large compared to their radii. The interest in a thermostatted electric field 
arises from the fact that at small fields a transport coefficient, the electrical conductivity of 
the particles, is proportional to the sum of the Lyapunov exponents describing the chaotic 
motion of the moving particle The Lyapunov exponents are to be calculated for the 
case where the charged particle is described by a non-equilibrium steady state phase-space 
distribution function which is reached from some typical initial distribution function after 
a sufficiently long period of time. In this state, the distribution function for an ensemble 
of moving particles (all interacting with the scatterers and the field, but not with each 
other) is independent of time and its average over the distribution of scatterers is spatially 
homogeneous. It is known from computer simulations [3] and theoretical discussions ||. ^[ 
that in the stationary state the trajectories of the moving particles in phase space lie on a 
fractal attractor of lower dimension than the dimension of the constant energy surface, which 
is three dimensional for the constant energy Lorentz gas in two dimensions. There can be 
at most two non-zero Lyapunov exponents for this model since the Lyapunov exponent in 
the direction of the phase-space trajectory is zero. Also, the relation between the Lyapunov 
exponents and the electrical conductivity requires that the sum of the non-zero exponents 
should be negative due to the positivity of electrical conductivity 

The case of the dilute, random Lorentz gas has already been studied in detail. Van 
Beijeren and coworkers [7| have calculated the Lyapunov spectrum for an equilibrium 
Lorentz gas in two and three dimensions using various kinetic theory methods including 
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Boltzmann equation techniques. These methods were also applied to the dilute, random 
Lorentz gas in a thermostatted electric field with results for the Lyapunov exponents that 
are in excellent agreement with computer simulations || ||. Moreover, the results for the 
field-dependent case were in accord with the relation between the electrical conductivity and 
the Lyapunov exponents for the moving particle. 

The purpose of this paper is to extend the results obtained for the Lyapunov exponents 
for dilute Lorentz gases to higher densities. Our central themes will be: (a) to describe a 
general method, based upon the BBGKY hierarchy equations, for accomplishing this task, 
and (b) to examine the effects on the Lyapunov exponents of long range in time correlations 
between the moving particle and the scatterers produced by correlated collision sequences 
where the particle collides with a given scatterer more than once and the time interval 
between such re-collisions is on the order of several mean free times, with an arbitrary 
number of intermediate collisions with other scatterers. These correlated collision sequences 
are of particular interest in kinetic and transport theory because they are responsible for the 
"long-time-tail" effects in the Green-Kubo time correlation functions, which lead to various 
divergences in the transport coefficients for two and three dimensional gases, where all of the 



particles move [10|. In the case of a Lorentz gas in d dimensions, the Green-Kubo velocity 



correlation functions decay with time, t, as t - ( d / 2+1 ) [11] and the diffusion coefficient is finite 
in both two and three dimensions. Here we describe the effects of these type of correlations on 
the Lyapunov exponents for the two-dimensional Lorentz gas, in equilibrium, where we find 
no effect, and in a thermostatted electric field, where we find a small, logarithmic dependence 
of the Lyapunov exponents upon the applied field. This logarithmic effect is an indicator for 
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similar effects to be expected when one calculates Lyapunov exponents associated with more 
general transport in two-dimensional gases, whereas in three dimensional systems one would 
expect corresponding non-analytic terms proportional to e 5 / 2 . In the case of the Lorentz 
gas, at least, the logarithmic terms can easily be associated with the logarithmic terms that 
appear in the field dependent collision frequency, and a very simple argument can be used 
to establish this relation between logarithmic terms in the Lyapunov exponents and in the 
collision frequency 

In Section 2 of this article we describe the general theory of Lyapunov exponents of a two- 
dimensional thermostatted electric field-driven Lorentz gas and quote the results within the 
scope of the Boltzmann equation. In Section 3, we generalize the theory to incorporate the 
effect of correlated collision sequences. In Section 4, we outline the calculation of the effects 
of the correlated collision sequences on the non-zero Lyapunov exponents using the BBGKY 
equations discussed in Section 3 and obtain the non-analytic field-dependent term in the 
Lyapunov exponents, originating from the correlated collision sequences, along with other 
analytic field-dependent terms. In Section 5, we present some simple arguments explaining 
the field- dependence of the collision frequency and show that this is the sole origin of the non- 
analytic, field-dependent terms in the Lyapunov exponents. Notice that the arguments given 
in Section 5 are independent of and much simpler to follow than the formalism developed in 
Sections 3 and 4. We conclude in Section 6 with a discussion of the results obtained here, and 
with a consideration of open questions. Methods for determination of the field-dependence 
of the collision frequency are outlined in the Appendix. 
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2 Lyapunov exponents of field driven Lorentz gases in 
two dimensions 

2.1 General theory 

The random Lorentz gas consists of point particles of mass m and charge q moving in a 
random array of fixed scatterers. In two dimensions, each scatterer is a hard disk of radius 
a. The disks do not overlap with each other and are distributed with number density n, 
such that at low density na 2 < 1. The point particles are acted upon by a uniform, constant 
electric field E in the x direction, but there is no interaction between any two point particles. 
There is also a Gaussian thermostat in the system to keep the speed of each particle constant 
at v by means of a dynamical friction during flights between collisions with the scatterers. 
The collisions between a point particle and the scatterers are instantaneous, specular and 
elastic. During a flight, the equations of motion of a point particle are 



r = 



P 

v = — , 
m 



p = mv = qE — ap 





v_= (v cos9_,v sin0_) 



v + = (v cos0 + , v sin0 + ) 



Fig. 1 : Collision between a point particle and a scatterer. 
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and at a collision with a scatterer, the post-collisional position and velocity, r+ and v+, are 
related to the pre-collisional position and velocity, r_ and V-, by 

r + = r _ , tT + = iL — 2 • <r) a , (2) 

where a is the unit vector from the center of the scatterer to the point of collision (see Fig. 
1). The fact that each particle has a constant speed v determines the value of a : 

a = — => p = qE —p. (3) 

Equivalently, in polar coordinates, the velocity direction with respect to the field, defined 
through v • x — cos 9, changes between collisions as 

= -esinO, (4) 

where e = and we define the dimensionless electric field parameter e — — , where 

mv v 

I = (2na)~ 1 is the mean free path length for the particle in the dilute Lorentz gas. To 
denote the electric field, we will normally use e, though from time to time we will use e, too. 

Treating this two-dimensional Lorentz gas as a dynamical system, we define the Lyapunov 
exponents in the usual way: a point particle in its phase space (r, v) — X starts at time 
to at a phase space location X(t )- Under time evolution, X(t) follows a trajectory in 
this phase space which we call the "reference trajectory". We consider an infinitesimally 
displaced trajectory which starts at the same time t , but at X'(t ) = X(t ) + 5X(to)- 
Under time evolution, X'(t) follows another trajectory, always staying infinitesimally close 
to the reference trajectory. This trajectory we call the "adjacent trajectory". Typically the 
two trajectories will separate in time due to the convex nature of the collisions. Thus, we 



can define the positive Lyapunov exponent as 

A+= lim , im IJ^ + ^l, (5 ) 

T^oo |5X(t )H0 T |5X(t )| 

for a typical trajectory of the system. 

We assume that, for small fields, this Lorentz gas system is hyperbolic. Since the two- 
dimensional Lorentz gas can have at most two nonzero Lyapunov exponents, we denote 
the negative Lyapunov exponent by A_. Without any loss of generality, we can choose to 
measure the separation of the reference and adjacent trajectories equivalently in r-space, 
thereby reducing the definition of the positive Lyapunov exponent to 

1 \6r(t + T)| 

A+ = hm hm -In — — — — . (6) 

T^oo |5r(t )H0 1 I^Wl 

In order to calculate the right hand side of Eq. we introduce another dynamical quantity, 
the radius of curvature p, characterizing the spatial separation of the two trajectories (see 
Fig. 2) : 




In Fig. 2, a particle on the reference trajectory would be at point A at time t. At the 
same time, a particle on the adjacent trajectory would be at B. A local perpendicular on the 
reference trajectory at A intersects the adjacent trajectory at C. The backward extensions 
of instantaneous velocity directions on the reference and adjacent trajectories at A and C, 
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respectively, intersect each other at point P. We denote the length of the line segment AC 
by SS(t) and ZAPC by 59(t). The radius of curvature associated with the particle on the 
reference trajectory at time t is then given by 

„M = f§ = l^|. (7) 

Having defined p(t), one can make a simple geometric argument to show that 

8S(t) = v59{t) = (8) 

so as to obtain a version of Sinai's formula [[EJ, 

r to + T a + 
v at 



A = lim - / — -. (9) 

T^oo T J. p(t) 



to 

During a flight, the equation of motion for p is given by || 

p 2 e 2 sm 2 9 , . 

p = v + pecosO + . (10) 

v 

At a collision with a scatterer, the post-collisional velocity angle 9 + and radius of curvature 
p + are related to the pre-collisional velocity angle 9- and radius of curvature p_ by |13], [T3 



11 2 e 

9 + = 9- - 7T + 2(fi, — = 1 h - tan0(sin6'_ + sin^ + ) , (11) 

p + p- acoscp v 

where cj) is the collision angle, i.e, cos0 = \v- ■ a\ = \v + ■ a\ (see Fig. 1). 

Now we assume that, for sufficiently weak electric field, the field-driven Lorentz gas in 
two dimensions is ergodic, and that we can replace the long time average in Eq. (|B|) by a non- 
equilibrium steady state (NESS) average, including an average over all allowed configurations 
of scatterers, to obtain 



A=(-\ ■ (12) 



PI 



NESS 
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The electric field is considered weak if the work done by the electric field on the point particle 

over a flight of one mean free path is much smaller than the particle's kinetic energy, i.e, 
q\E\l el 

-'—5- = — =£<< 1. 

We note for future reference, that the sum of the two nonzero Lyapunov exponents is 
related to the average of the friction coefficient a, by [l], [15|, [16 



,'qE-v\ J-E oE 2 . . 

A + + A - = -<«>N ESS = -( 1 -r) =-^ = -— r ( 13 ) 



mV ' NESS mV mV 



Here the electric current J = (qv) NESS is, for small fields, assumed to satisfy Ohm's law, 
J = aE, and a is the electrical conductivity. 



2.2 Results obtained using the Lorentz-Boltzmann equation 

To the lowest order in density, one can assume that the collisions suffered by the point 
particle are uncorrelated, and use an extended Lorentz-Boltzmann equation (ELBE) for 
the distribution function of the moving particle, fi(f,v,p,t) in (r, v, p)-space [|J needed for 
the evaluation of the averages appearing in Eqs. ( |T2D and ([13]). To calculate the positive 
Lyapunov exponent, one needs to consider the forward-time ELBE while to calculate the 
negative Lyapunov exponent one needs the time reversed ELBE. To the leading order in 
density, the Lyapunov exponents are then given by f| : 

Af = Ac - §-e 2 + 0(e 4 ) and AL B) = - A - ^--e 2 + 0(e 4 ) . (14) 
48 v 48 v 

The superscript, B, indicates that these are results obtained from the Lorentz-Boltzmann 
equation. Here A is the positive Lyapunov exponent for a field-free Lorentz gas (see for 
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example |I7|]) given by 

X = 2nav[l-C -\n(2na 2 )}, (15) 



where C is Euler's constant, C = 0.5772.... From Eqs. (|T3|) and (|T4|), using Einstein's relation 



between diffusion constant and conductivity, one gets the correct diffusion coefficient within 

3 

the Boltzmann regime, D (B) = - Iv. 

8 

To derive the results in Eq. (|14D , one uses Eq. (|TT| ) with e = 0. The ^-dependent term in 
Eq. ( jTTD can be explicitly shown to be of higher order in the density than the terms present in 
Eq. (|14D [18] . In the following sections we will investigate the effect of sequences of correlated 



collisions between the point particle and the scatterers on the Lyapunov exponents. However, 
the e-dependent term in Eq. ( |TT]) will again be neglected since we will present the effect of 
these correlated collision sequences in leading order in the density of scatterers only. Thus, 
instead of Eq. flll]) , we will use 

112 

— = — + -. 16 

p + p_ a cos 

3 The extension of the ELBE to higher density 

3.1 Binary collision operators in (r, v, p)-space and 

the BBGKY hierarchy equations 

The Boltzmann theory for the Lyapunov exponents assumes that the scatterers form a dilute, 
but quenched system and that the collisions of the point particles with the scatterers are 
uncorrelated. To incorporate the effects of correlated collisions on the Lyapunov exponents, 
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we will use a method based on the BBGKY hierarchy equations, familiar from the kinetic 



theory of moderately dense gases [jnj. Since the moving particles do not interact with each 
other, it is sufficient to consider the distribution functions for just one of them, together 
with a number of scatterers. One starts from a fundamental equation for an (N + l)-body 
distribution function, /n+i = fN+i(r, v, p; Ri, R 2 , ., ., Rn', t), which is the probability density 
function in the entire extended phase space T spanned by the variables r, v, p, Ri, R 2 , ., ., Rn, 
describing our system of N scatterers and one moving particle. We require that /jv+i satisfies 
the normalization condition 

dfdvdpdRidR 2 --dRNfN + i(f,v,p;Ri,R2,.,.,R]y;t) = 1. (17) 

This (N + l)-body distribution function satisfies a Liouville-like equation determined by the 
collisions of the moving particles with the scatterers and by the motion of the particles in the 
thermostatted electric field, between collisions. Since the time evolution of r and v in this 
field is not Hamiltonian, we must use the Liouville equation in the form of a conservation 
law, rather than the usual form for Hamiltonian systems, to obtain 

+ %'(?fN + i) + V^-^/iv+i) + ^-(p/iv+i) = 53f_, < / JV+1 . (18) 

i = 1 

Here the operators JL^ are binary collision operators which describe the effects on the 
distribution function due to an instantaneous, elastic collision between the moving particle 
and the scatterer labeled by the index i. The explicit form of the binary collision operators 



may be easily obtained by a slight modification of the methods used by Ernst et al. |L9 



order to include the radius of curvature as an additional variable. One finds that the action 
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of this operator on any function /(r, v, p; R 2 , ., ., -R,; £) is 



T_ - f = a da; \v • &A \ / dp' S\ p — — - - C ° S — - — | 5(aai — if — Ri)) x 

x b aiipl - d^aa, + (f - £))}/, (19) 

where <jj is the unit vector from the center of the scatterer fixed at Ri to the point of collision. 
The action of the operator b auP > on the function f(r, p; R%, R2, ., ., -R?; t) is defined by 

ba i>P ' f(r,v,p;R u R 2 ,.,.,Rj;t) = f(f,v-2(v-a i )a il p';R h R 2l .,.,Rj;t); (20) 

that is, b aiP i is a substitution operator that replaces p by p' and the velocity v by its restituting 
value, i.e, the value it should have before collision so as to lead to the value v after collision. 
It is often useful to express the binary collision operators as as sum of two terms such that 

T-,i = f^l - TH , (21) 

where 

= a [ ddi\v-ai\ [ dp'sip - , C ° S ^' - ] 5(aai - (r - Ri))b ai:p/ (22) 

i^>o Jo V « cos <^ + V/ 

and 

fi'j = a I ddi \v -bi\ 5{aai + (r - k\)) . (23) 

J V-(7i > 

One sees that T^J f and t1~^ / respectively describe the rate of "gain" and the rate of "loss" 
of / due to a collision of the point particle with the scatterer fixed at Ri. 

The BBGKY hierarchy equations are then obtained from Eq. ([18]) by integrating over 
scatterer coordinates, as a set of equations for the reduced distributions for the moving 
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particle and (j — 1) scatterers, denned by 



N\ 



-yyy J (IRj . AR N f N + i (f , v, p; R t , R 2 , . , . , R N ] t) . (24) 



(N-j + l)\ 

One then easily obtains the BBGKY hierarchy equations (1 < j < N) 



at 



(25) 



3.2 Cluster expansions and truncation of the hierarchy equations 

The usual procedure for truncating the hierarchy equations in order to obtain the Boltzmann 
equation and its extension to higher densities is to make cluster expansions of the distribution 
functions, f 2 , ■ ■ ■ in terms of a set of correlation functions, g 2 , g% ■ ■ ■ as follows: 

f 2 (r,v 1 p;R 1 ;t) = nfi(f,v,p;t) + g 2 (r,€,p;Ri;t), (26) 

hif,v,p]R u R 2 ]t) = n 2 fi(r,v,p]t) + ng 2 {r,v, p;R x ;t) + ng 2 (r,v, p; R 2 ;t) 

+ g^r^p-RtiR^t), (27) 

and so on. Hereafter, to save writing, we denote g 2 (f,v,p;Ri;t) as g 2 ^ , g 2 (r,v, p; R 2 ;t) 
as g 2 /s(r, v, p; R\,R 2 ; t) as fs and g% (r, v, p; Ri,R 2 ; t) as #3. The first terms in each of 
these expansions represent the totally uncorrelated situation, where there are independent 
probabilities of finding the moving particle and the scatterers at the designated coordinates. 
The next terms involving the pair correlation functions g 2 ^ in Eqs. (|26| ) and (|2T| ) take into 
account possible dynamical and excluded volume correlations between the point particle and 
the scatterer at If one replaces f 2 by nfi in the first BBGKY hierarchy equation, Eq. 
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with j = 1, reduces to the ELBE. To find the corrections to the ELBE for higher densities, 
one must keep the g 2 term in Eq. (fffi ) and use the second hierarchy equation to determine g 2 - 
However, in order to solve the second equation, we have to say something about g^. A careful 
examination of the second and higher equations shows that g% contains, of course, the effects 
of three-body correlations, i.e, correlated collisions involving the point particle, a scatterer 
fixed at Ri and another scatterer fixed at R 2 , as well as excluded volume corrections due to 
the non-overlapping property of the scatterers. Here we will be primarily interested in the 
effects of the so called "ring" collisions on the Lyapunov exponents. These collision sequences 
are composed of one collision of the moving particle with a given scatterer, followed by an 
arbitrary number of collisions with a succession of different scatterers, and completed by a 
final re-collision of the moving particle with the first scatterer in the sequence, as illustrated 
in Fig 3. 




Fig. 3 : Sequential collisions with scatterers at R2, R3, ., . adding up to the ring diagram. 



The ring diagrams, taken individually, are the most divergent terms that appear in the 
expansion of dynamical properties of the Lorentz gas as a power series in the density of 
scatterers. They lead to the logarithmic terms in the density expansion of the diffusion 
coefficient of the moving particle |50|, and to the algebraic long time tails in the velocity 
time correlation function of the moving particle [jnfl . While many other dynamical events 
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and excluded volume effects contribute to the Lyapunov exponents, and must be included 
for a full treatment, we concentrate here on the effects of these most divergent collision 
sequences, since in other contexts, they are responsible for the most dramatic higher density 
corrections to the Boltzmann equation results. 

Thus we drop g% in Eq. (|27f) and obtain a somewhat simplified cluster expansion of 
given by 

h = n 2 h + ng 2 Ri + ng 2 ^ (28) 

Using Eqs. (|6|) and ( [28D and the first two of the BBGKY hierarchy equations, we obtain a 
closed set of two equations involving two unknowns, fi and g 2 $ given by 

V* • (CT ft) + — {ph) = J dEif_ fl [nh + <? M J (29) 

and 

- d f - ~ 

(^2,^) + v ^>- (^2,^) + -Q-piPQ^nd ~ n J dR 2 T_, 2 g 2 ni = nT_ (30) 

In the derivation of Eq. (^) from the second hierarchy equation not only have we dropped 
#3 as discussed above, we also dropped a term of the form T_ ig 2 g . This term provides 
"repeated ring" corrections to the ring contributions to g 2 g i . These are of the same order 



as terms neglected by dropping g 3 [^T], [L8[ , and should be neglected for consistency. We also 
dropped the time derivatives in the equations, so we are now looking for the distribution 
and correlation functions appropriate for the NESS. 

In Section 4, we will solve Eqs. ( |2DD and ( ppp in order to calculate the ring contributions 
to the positive Lyapunov exponent. Before doing so, it is useful to write down the usual 
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form of the ring equations in (r, -y)-space, which can be obtained by integrating Eqs. ( p9| ) 
and (POD over all values of the radius of curvature, < p < oo. We define the usual 
single-particle distribution function by, F\ = J p>Q dp f% and the pair-correlation function 
^2 = f p >o dp 92 r x ■ By imposing the boundary conditions that both fi and g 2 ^ go to 
zero as p — > and as p — > oo, we obtain 



V*-(tfiq) = y 4t_,i [nFi + G 2 /? J 



(31) 



and 



Vf(fG ai ) + V^-(#G 2j ^) - nj dR 2 T_ !2 G 2jii = nT_ A F x . (32) 

The actions of T"_ i or T_ 2 on Fi and G 2 ^ can be obtained by appropriately integrating 
T_ i /i, T-,1 <? 2 ^ or T_ )2 # 2 $ over p from to oo using the definitions in Eqs. ( |TP| ) and 
(P0|). T_ ( i and T_ 2 are the analogs in (r, -y) space of T_ ! and T_ 2 (see Eqs. ([ID]) and 
©), i.e., 

TL, = a / d<7j |u • o-j| < 5 (aa-» - (f - ^)) 6 CTi - 5(a&i + (f - fly) ^ . (33) 

In future applications however, we will drop the adi terms from the arguments of both 
5 (a&i ± (r — in T_ 5 j and T_ j: ; operators since they lead to corrections similar to 
excluded volume terms, neglected already. 



4 Effects of long range time correlation on A + and A_ 

We now concentrate on the solution of the BBGKY equations for the distribution functions 
that determine the Lyapunov exponents. The solutions of Eqs. (|2D| ) and (|3T)| ) are to be 
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obtained as expansions in two small variables, na 2 and e. The density expansion will give 
the corrections to the previously obtained Boltzmann regime results from the ELBE, and 
the e expansion will provide the field dependence of these corrections. We therefore write 
the density expansions of fx and g 2 (hereafter we drop the subscript Ri from g 2 ^) as 

fx = ti B) + ti R) + ■■■ and g 2 = + . . . , (34) 

where the superscript B indicates the lowest density result for fx as given by the ELBE, and 
the superscript R denotes the ring contribution. At the order in density of interest here, the 
quantities indicated explicitly in the above equations satisfy 

V*- (iT/D + §- p (pfi B) ) = nj dRxf^fP, (35) 

v*- (tfn + §- p (pfn = J dR.f^infr + <?n (36) 

and 

V r --(r^ R) ) + %'&gF) + ^P9? ] ) ~nj dR 2 f^ 2 g™ = nf^f™. (37) 

Our aim here is to solve Eqs. (|36| ) and fl37|) using the results of the ELBE for . We 
suppose further that each of these functions possesses an expansion in powers of e as 

A (B ' R) = AT + eflT + ^/i (B 2 R) + ••• (38) 

and 

9?" = 9^ + eg^ + e 2 g^ 2 + ... (39) 

The functions f^l have been previously obtained as the e solutions of the ELBE. Since we 
will be dealing with g 2 only in the context of of the ring term, we drop the superscript R 
from now on. 
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As mentioned above, we will neglect the term aa within the arguments of the 5-functions 
appearing in each of the binary collision operators T_ and T_, so as to take the moving 
particle to be located at the same point as the center of the appropriate scatterer at collision. 
The terms neglected by this approximation lead to higher density corrections to the terms 
we will obtain below. Secondly, an inspection of the radius of curvature delta function in 
the expression for the "gain" part of the binary collision operator, Eq. (p2|), shows that 

this term is only no n- vanishing when p < — , and that the dominant contribution to the p' 

r ! / i at n p' a cos (pi a cos (pi, ^ . 

integration comes from the region p ~ I. Naturally, ■ ~ (1 + Oynj) 

a cos (pi + 2p' 2 

in the argument of the delta function. In the Boltzmann level approximation this 0(n) 
term may therefore be neglected and it can be shown not to contribute to the leading field- 
dependent ring term effects on the Lyapunov exponents. Therefore, we will neglect it in 
what follows. Under this approximation, the gain part of the binary collision operator || 
acts on an arbitrary function h{f, v, p) as 

f^ h(f,v,p) « 5{r- Ri) - p)I(p)[H(r,v + ) + H(r,vj] = S(f— Ri) T(p, H) .(40) 
Here 

v± = v - 2 (v ■ a it± ) & ij± , (41) 
and &i t ± is defined by the condition that the scattering angle (p = ±cos _1 (— j. Also 



oo 

H(r,v) = I dp'h(r,v,p') (42) 
Jo 

and 

,(„) = , 4 "" , . (43) 
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Finally, we express the velocity vector v in terms of the angle that it makes with the 
direction of the electric field so as to obtain the following set of equations for the terms in 
the ^-expansion of g 2 



v ■ Vf + v— h 2nav 

Op 



g 2 , = -2nav5(f- R x ) f™ + n5{?-R x )Y{p, F™) , (44) 



v ■ Vf + v— h 2nav 

Op 



9 ( „ \ d 
92,i = ~ ^Hp cos 9g 2 ,o I + — I sin # g 2> o 

-2nau<y(r - #i) + n 5(r - R x ) T(p, F$) , (45) 



and 



v • + v— — h 2nai> 
ap 



9 / „ p 2 sin 2 6> \ d 
92,2 = - q- ( P cos 6/ g 2 , i H #2,0 ) + ^ I SUI 1 



- 2nav S(r - R\) f[ B) 2 + n 5{r- R\) T(p, F^) . (46) 



We notice that the equations thus obtained are linear inhomogeneous differential equations 

_ Q 

of the form L g 2 j = bj (j = 0, 1, 2), where L = Vf + v— — h 2nav~\ is a linear differential 
operator and bj for j = 0, 1, 2 are the inhomogeneous terms on the r.h.s. of Eqs. (fHj), 
and fl46|) respectively 



We will need to solve Eqs. (§4|-j46|), in conjunction with the equations for G 2i o, G 2} i and 
G 2}2 obtained by directly integrating Eqs. ([44l-|4"6"|) over p from to 00. The equations for 
the corresponding G 2 $ are then 



v-Vjt—n I dRo T 



1-2 j- -,2 



G*2,o — nT^^i F-^ , 



(47) 



v ■ Vf — n dR 2 T 



G 2 , 1 



d_ 

m 



sm9G 2 , ) + nT^F™ 



(4£ 
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and 



v-V? - nj ' dRnT.,2 G 2 , 2 = fsin0G 2)1 ) + nT. tl F$ . (49) 

The equations for G 2} q, ^2,1 and ^2,2 are also linear differential equations of the form 
L G 2 j = Bj, where l! — [t? '• Vf — n J a?i?2 ^-,2] and -Bj- for j = 0, 1, 2 are the inhomoge- 



neous terms on the r.h.s. of Eqs. p7|) , (|4q) and (|49|) respectively. 

To solve these equations, we take Fourier transforms of g 2 and of G 2 in the variables r and 

in the velocity angle, 9. That is, we define g 2 (k) = — = [ dr g 2 e~ ~ ^ and calculate 

yV Jv 

92,o{k), g2,i(k) and g 2 , 2 (k) in the fc-basis, using periodic boundary conditions. Similarly, we 
define the m-th angular mode of g 2 (k) as g 2 m \k) = / d6 e~ ~ %m6 g 2 (k) . We also define 

V27T J 

G 2 m \k) in an analogous way. Thus, corresponding to Eqs. ( f44| - f46D and Eqs. ([47)-|4"9"|) , we 
have two sets of three equations to be solved, one involving g 2> o(k), g2,i(k) and g 2 , 2 (k), and 
the other involving G 2> o(k), G 2} i(k) and G 2j2 (k), in (k, m) basis. In this basis, the operators 

Q _ _ 

Lj: = [ik-v + v— — h 2navj and = [ik-v — n f dR 2 T^ t2 ] are both infinite dimensional 
matrices in m-space and both of them have non-zero off-diagonal elements due to the term 
ik ■ v generated from the operator v ■ Vf . However, it is easily seen that these off-diagonal 
elements are proportional to <5 m , m +i and <5 mim _i and they are easily treated. 

A further simplification can be made by noticing that the schematic forms of the solutions 
are g 2 j(k) = [L^ 1 bj(k) and G 2 j(k) = [ifj _1 Bj(k) and hence the dominant parts of g 2 ,j(k) 
and G 2 j(k) will come, loosely speaking, from the eigenfunctions of and L-. having the 



smallest eigenvalues. The lowest eigenvalues of L t are oc k 2 due to the contributions from 



the hydrodynamic modes |TI|. Thus, to capture the dominant part of the solutions we 



should solve the equations in the range k = \k\ « I x , the inverse mean free path and use 
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perturbation expansions in the small parameter kl. We will not, in our analysis, follow the 
mode expansion technique, as it is simpler to calculate G 2 directly. However, one can use 
mode expansions and one finds that the results of both the methods agree. 

4.1 Solution for G2 

To solve for G 2 first we need to know the solutions of the Lorentz-Boltzmann equation for 
F^q, Fl B l and F$. These are given by § 

F$ = F$ = — ^— cos 6 and F$ = f cos 29. (50) 

We note that in m-space, defined above, the m-th diagonal element of the infinite matrix 

. 4:171^ V 

L r is — — -- while the off-diagonal elements are ikv 5 mm ±i- Thus an expansion in 

k [Am 1 — 1) / 

k = kl can be easily obtained by considering successively larger parts of the matrix L'^ in the 
index m, starting with 3x3, 5x5 matrices and so on, chosen in such a way that the element 
of corresponding to m = appears as the center element of these matrices. As we want 
to make our results correct up to 0(k°), we need to increase the size of these matrices till the 
expressions of G 2 j{k) obtained from G 2 j(k) = [-^r] -1 Bj(k) (f° r j — 0, 1, 2) converges up to 
0(k ). Also, as we want to obtain the expression of A + and A_ in the leading field-dependent 
order, which is e 2 ; we need the solutions of all the m-modes of G , 2,o(^), G2,\{k) and (j2,2(^) 
that are necessary to obtain all the terms of /{ R) that are oc e 2 and contribute to this leading 
field-dependent order of A + and A_. In more explicit form, this means that we definitely 
need the solutions of G { 2 m = °\k), G { 2 m = \k), G { 2 m = ±1] (k), G ( 2 m = ±2 \k) and G ( 2 m 2 =0) (k) up to 
0(k°). However, once we present these solutions, from the structure and properties of them, 
it will turn out that we will also need the expressions of G 2 m 2 ~°\k) in the leading order of k 
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for j = 2, 3, . ., to consistently obtain all the terms, that are oc e 2 . 
At the e° or equilibrium order, we find 

G { $(k) = Vm. 

Proceeding to order e, we find that Gr^i (A;)'s obey 



(51) 



IV 



(k X + lk y )G { ™+ l) {k) + (k X - lky)G^-'\k) 



(m-1). 



+ 



inavm- ^ (m) 



4m 2 



2 V2Ty 



{8 m ,i + 5m, -l) ,(52) 



where it turns out that we need a 5 x 5 matrix block corresponding to m = —2, —1, 0, 1, 
and 2 to get the solutions of G\ ™ (k) up to 0(k°) that are relevant for us, yielding 



Gtr°\k) 



ctr l \k) = - 



IQnavk 2 y/2irV 



Gtr 2 \k) 



45 ky k% iky ^ 



Gtr l \k) 



and G^~ 2 \k) 



lQnavk 2 ^/2nV 



45 ky [k>x ^ky 



(53) 



1024(na) 2 vk 2 ^2W ' 1024(na) 2 vk 2 V2nV 

Notice that we have also calculated G^i~ 2 \k) and &£]~ 2 \k) } even though they are O(k), 
because they affect the 0(k°) solution for G^^ \k)- 

For order e 2 , the relevant G^ (^)' s are then calculated using Eq. ( |5~3"D and considering a 
5x5 matrix block of L'-.. There we need only the solution for G^™ -0 ^) : 

kl 45 (2k 2 x - 5k 2 ) 



r 2,2 



v 2 A; 4 v / 2W 1024(rMw) 2 £; 2 



(54) 



Examining the properties of the solutions, Eqs. ( [53] - |54] ) and observing from Eqs. (^7|44D|) the 
way the solution of G2J affects the solution of G 2 , (j + i), one sees that the leading power of k 
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in the expression of G^^°\k) for (j = 1,2,3....) is k~ 2 K However, in the expression of G 2 , 
&™2^°\k) appears with a factor of e 2 K When G2 is finally calculated, after a summation 
of the appropriate &-values[], the contribution of the sum of all the effects coming from the 
Q(k- 2 i) terms of the G^J 0) (jb)'s is seen to be in the same order of density of scatterers as 
the 0(k°) term on the r.h.s. of Eq. (p)4l). In fact, it also turns out that the 0(k~ 2 i) terms of 
the G 2 w ^ Q \k)''s are the only ones among the G 2 m j(kys that contribute to f[ R) in the order 
of e 2 . This implies that along with the solutions, Eqs. (|5lD , (|53|) and (|54D , we also need to 
include the 0(k~ 2: >) term of G 2 m 2 j°\k) to be consistent. If one just considers this 0(k~ 2j ) 

V(m=0), 



term in G 2 2 j then it is easy to see that they satisfy a recurrence relation for j > 1 

■-i(m=0) /A _ r«x Mm= 

r 2,2(i + l)l ft J - v 2 k 4 ^2,2j 



~( 0) " k 2 ~r 0) 



2 ' 2j 1 ; V^ 4 



The solutions, Eqs. (pTj), (|53"D, (0) an< ^ are then used to determine the integration 
constants that arise when we solve the differential Eqs. fl44l - f46|) . It is important to note that 



the first term on the right hand side of Eq. fl54|) is inversely proportional to k 2 . This is the 
origin of the logarithmic terms we find below. 



4.2 Solution for g 2 

Here we apply the same procedure to solve for the #2 (fc)'s from the equations = bj(k) 

for j = 0, 1, 2. This time, the elements of are differential operators in the variable p and 
1 To see how the fc-integration is performed, see the last paragraph of Section 4.2 



24 



the corresponding constants of integrations are determined using the solutions of (5 2 m ^(fc)'s 
while maintaining that g% (&)'s go to zero as p — > and as p —>■ 00. We also note that 
for our purpose, solutions of the city's are only needed for p > — as the solution of the 
(^(^O's for p < ^ gives rise to higher order density corrections than under consideration here. 
These solutions can also be obtained by the mode expansion technique discussed above in 
the paragraph preceding Section 4.1. However, as it is fairly straightforward to solve the 



differential Eqs. (44-46) for p > — , we directly write down the necessary solutions up to 



0(k°). 



We obtain, for p > 



2na 
V2W 



2nap) e 



- 2nap 



_( m =-l) 



(k) 



k x {k x ~\~ iky) ^ k x {k x + iky) 2 n< ^p _|_ ^ na, P 



8k 2 



2k 2 



+ 



(2nap) 2 (2napY 



■ Inap 



~9ir°\k) 



)k 2 V2%V 



2nae~ 2nap , 



#i =1 ¥) 



■y/2nV 



k x {k x iky) k x (k x iky) 



8k 2 



+ 



2k 2 



2nap + 



2nap 



+ 



(2nap) 2 (2napf 



- 2nap 



v 2 k*V2W 
1 



2nae' 2nap 
45 



2v 2 V2nV 



315^ hi 11 2 

"77 + 777 P - 77 2na P 



128na 256nak 2 Ak 2 

13 ^ 2 19, , 2 „3 1 /r.„\2„3 



8 /c 2 



272a// + — (2na) Y + — ^ (2na) > 
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2/c 2 



1 "\ 1 
+ ±£ (2na)V - - (2na) 4 p 5 
24 v ; H 8 V ; P 



-2nap 



(57) 
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and for the 0(k 2j ) terms in (k) we have 



gi m =°\k) = { -4^(4l] J ^ae-^. (58) 



We point out that all of the terms in each of the square brackets, in each of the above three 
equations, are of the same order in the density. This can be seen easily by noting that p 
is typically of order (2na) _1 , so that (2nap) is typically independent of the density. The 
solutions, Eqs. (|5l|), (|53D , (|54|), (|56|), d57|) and ([58]) now can be assembled to calculate G% 



and g 2 in (r, v) and (f,v,p) space respectively and feed the results into the r.h.s. of Eqs. 
(|3T|) and fl36|) to obtain f[ R> . This involves a summation of different m and A;- values. In the 



infinite volume limit the /c-sum can be converted to an integration over k. The sum over m 

is straightforward, but we have to remember that the integration over k has to be carried 

out in a range k < k ~ l^ 1 . Secondly, since we have expanded the distribution functions in 

powers of e and then subsequently in powers of k, the lower limit of k for the ^-integration 

cannot be taken to be zero. To determine this lower limit of k for the /c-integration, we 

observe that the expansion in e cannot be carried out for those values of k where k < — , so 

that the value — forms a natural lower cut-off for the Fourier transform. Our solutions of 
2v 

G2 and #2 therefore do not hold for k < — and to do a satisfactory perturbation theory in 

2v 

the range k < — , one needs to consider both the e and ^-dependent terms together. After 
2v 

doing so, one finds that such a perturbation theory does not affect our results at the present 



density order [18 



Before performing the integration over k, we notice that in two dimensions, the numerator 
of the fc-integral is proportional to k dk. This means that any part of the solutions of 
cj2{k) or G2{k) having a leading power of k of order 2 or higher in the denominator gives 
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rise to a singularity at k — > for the k- integral. First, the highest leading power of k 
in the denominators of Eqs. (51), (f>3|), (54) and (p7|) is A; 2 , occurring in 

^T 0) (fc) and 



g^l respectively. These terms proportional to A; 2 give rise to a logarithmic electric 

-, s 
field dependence once the ^-integration is performed for — < k < k$. The rest of the terms 

2v 

in these solutions supply only analytic field dependences that can be expressed as power series 
in e. Secondly, even though the solutions given in Eqs. ( |56|) and (|58|) have higher powers of 
k than k 2 in the denominators, they also come with subsequently higher powers of e in their 
numerators. Thus, when the ^-integration is performed, they contribute terms proportional 
to e 2 or higher, to g 2 or G 2 - Consequently, in addition to analytic field dependent terms, in 
our present approximation we have only one non-analytic field dependent term appearing in 
#2 or G 2 and that is proportional to e 2 In e. No doubt there exist further non-analytic terms 
in higher orders in e, but their calculation would require a careful consideration of various 
terms we have neglected here, such as the repeated ring contributions. 

4.3 Solution for /{ R) and the calculation of A+ } 

Once the solutions, Eqs. (§TJ), Q53|), (0), ([56]), Q57D and fl5|) are inserted in Eqs. (0) and 
(|37j) and the ^-integration is performed in the range — < k < ko, we get, by the method 
described in (p]|-|4"5|) , the following equations to be solved to obtain /| R) and _F{ respectively 
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- e A (sin <,/<»>) + JLiL + p£cos0 + f <m\ + 2navf , 

= < ~ P) 



4:VP 

1 " (f) 1 



X 



— 2at> 



cos o CT 



in 2 v 2 



In 



2f A;q 



i-2 



8n 2 1 16nai> 



ecos6> 



135 



20A8(nav) 2 



2 1 



16ti 3 v 2 



in 2 v 2 



+ il e -2na P ( 2na (i _ 2nap) - - 
8n 2 v 



(2napf . N o 3 1 
_ (2nap) 2 - -2nap + - 



cos 9 



Anav 2 



(2nap) 1 



13 fn , 4 25 , n ,o 35 /ri , 2 2nap 135] 1 

- ^< 2 -rf 4 - M {2naf,) + i6 - -T + 25eJ) 



+ 



1 



16n 3 v 2 



As 2 2nae~ 2nap 



(59) 



and 



d 



— £ — (sin^Ff ') — nav / dcp cos (j) (b a — I) F{ 



09 



->(R) 



d0 cos (j) (b a — 1) 



16n 2 na 



as" . f'2vk 
In 



aw 



fcgf 3 



87r 2 \ 16naw 



ecosf? + 



135 



20A8(nav) 2 



167T 3 1 



As 2 + 



7^— e cos 9 + .... 



(60) 



where b a has been defined in Eq. (P3|). The ^-dependent terms in Eqs. (|59|) and (^0| 



originate from Eqs. ( |58D and flog ) respectively after the ^-integration is carried out. Here A 



is the integral^ 



.4 



2- 



cos 2 </> In 



1 + - cos 2 
4 



0.53536.... 



(61) 



2 We thank the referee for pointing out an error in a previous calculation of this integral. 
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The dominant effect of the ring term on the single particle distribution function, i.e, f[ 



(R) 

J 



can now be determined from Eqs. (59) and fl60|). It is also of some interest to give a crude 



estimate of the terms that we have neglected. One knows from other studies in the kinetic 



theory of gases |K| g2f that excluded volume corrections to Boltzmann equation results are 
the numerically most important corrections, until the density of the system becomes high 
enough that the mean free path of a particle is less than the size of the particle itself. These 
excluded volume corrections are provided by the Enskog theory, and this theory can be 



applied to the Lorentz gas, as well pOfl . In our case, the Enskog corrections can be included 
by replacing the density parameter n by n (1 — -nna 2 )^ 1 ~ n (1 + ima 2 ) in the Boltzmann 
equation. The Enskog correction affects both Ao and the e-dependent terms in the expressions 
for A± in Eq. fll^l). Along with the Enskog correction there are other correction terms that 



affect both A and the field-dependent terms A± JT8], Q at the same density order as the 
Enskog correction. Also, the terms that have been dropped to obtain Eq. (|TBp from Eq. 
(|TT|), contribute to A± at the same density order as the Enskog correction. However, since 
the principal objective of this paper is to investigate the non-analytic contribution of the 
ring term to the Lyapunov exponents, we will ignore the Enskog and related corrections from 
our consideration. Thus, using Eqs. ( p9[) and floTf), one can express the full solutions of fi 
and Fi as sums of a solution in the Boltzmann regime, a correction due to the ring term and 
a correction due to the Enskog term, plus all of the other terms we have neglected, as 

h = /r + /r + • ■ ■ ■ (62) 

F 1 = F/ B) + F™ + .... (63) 
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Consequently, for the positive Lyapunov exponent A + we have, 

A + = Af + Af + • • • • 
The solution of F X (R) is quite straightforward, 

F (R) 



(64) 



■ecosB + 



(65) 



128vr 2 (no) 2 f; 

However, to solve for we find that in addition to the analytic field- dependent terms 
which can be expressed as a power series in e, there is a non-analytic field-dependent term 
in /{ R) proportional to e 2 In e. Thus, with 



J 1 J 1 , analytic ' 



J 1, non-analytic 5 



(66) 



we have 

t/ 1, analytic 



2nap — 



(2napY 



+ 



4(nat>) i 



2nav 


I 4 


(2napf 


11 

+ 48 


32 



- 



79 



96 



3 . „ >o 135 , , 135 

+ 32 (2 " a " ) " 5l2 f 2 """' + 512 



+ 



(27r) 3 w 2 
135a/cg 
_(27r) 3 w 2 _ 512 (4™^ 



ry^A{l - 2nap)e- 2nap + . 



aA 



—2nap 

, a 
for p > - 

for p < — (67) 



and 



•/ 1, non-analytic 



47T 2 V 2 



47T 2 f 2 



In 
In 



2vfe 



2i>A; 



^1 — 2nap) e 



-2nap 



for p > 



s 



(?£)> 

v a 2 



for p < - . (6£ 



Notice that the ring contribution to the distribution function in Eqs. ( |6T| ) and satisfies 
the boundary conditions that /{ R) — > as p — > and p — >■ oo. Equations ( |5"7| ) and (|5"8| ) also 
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satisfy continuity at p — — at the leading density order. The distribution functions, Eqs. 
(|65|), (|37|) and (|68D, are all the ones that we need to calculate Consequently, 

WR) _ \ (R) , x(R) / fiQ \ 

/x + /V + , analytic ~T~ ^+ , non-analytic \ UJ / 

and using the the definition of Lyapunov exponents in Eq. (|12|) , we have 

MR 

(R) _ / ,/j / , Jl, analytic 



"27T (-00 f(B-) 

A^'analyt, = \ d9 d R 

J? P 



afcgu afc 2 /V f 13 135 



4?r 2tto 96 512 



(ln2na 2 + C) 



CLE^ 

0.53536- — ^- (In 2na 2 + C) + .... (70) 

(271")^ 



and 



/•2ir />oo 

A^non-analytic = d6 d P ^^ 

Jo h p 



|l{ ln ^} (ln W + C ), (71) 



where I is the mean free path and C is Euler's constant, C = 0.5772 .... 
4.4 Calculation of A w 

To calculate the corresponding effect of the ring term on A_, we make use of the relation 
Eq. flT3|). It is easy to calculate the effect of the ring term on («) NESS using already 
determined in the previous section. Thus, using 

X + + X - = - <«Lss > ( 72 ) 
and a complete analogy to Eqs. (|62] - |64D , we can calculate three terms of (tt) NESS : 



)nESS = < a >N B ESS + < a >N R ESS + ' " " " > ( 73 ) 
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with 



[a) m = 2^°L± + .... (74) 

v /ness 3 2vrw v 1 



Following Eqs. ( p2] - |64]) , we now express A_ as A_ = A_ + X^ + 

A (R) +A (R) /.\(R) 



-(«) 

\ It 



This leads us to 



satisfying 



\(R) 

— , analytic 



4:71 



2ttv 



5 

96 



135 
512 



ln2na 2 + C) 



+ 0.53536 



ae 



{2ir) 2 v 



In 2na 2 + C ) 



(75) 



and 



n_ lytic = ^ { In ^ } ( In 2na 2 + C ) . (76) 
where / is the mean free path and C is Euler's constant, C = 0.5772 .... 



5 The field-dependent collision frequency and its 
effects on the Lyapunov exponents 

As stated before, our second main purpose was the derivation of the leading non-analyticity in 
the field dependence of the Lyapunov exponents. In analogy with the transport coefficients, 
we expected these non-analyticities to result from the long time behavior of the ring terms, 
which we found confirmed in the preceding section. Some further thought reveals we can 
estimate the non-analytic field dependence in a simple way. 

In the presence of a thermostatted field there are two types of contributions to the positive 
Lyapunov exponent of the two-dimensional Lorentz gas: 
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1) contributions from the bending of the trajectories by the fields and 

2) contributions from the divergence of trajectory pairs at collisions. 

The first type of contributions are of order e 2 in the Boltzmann approximation. We 

expect that the coefficient of this term will pick up higher density corrections and there 

will be additional terms of higher orders in e. But we have not found any indications for 

corrections of lower order than e 2 resulting from the field-bending contributions. 

The collisional contributions can be generally expressed as an average of the form 
\Sv'\ 

v ( In — — - ) c , with Sv' and 5v the velocity differences between the adjacent trajectories just 
\bv\ 

after and just before a collision, respectively, v the average collision frequency, and the an- 
gular brackets, ( ) c , indicating an average over collisions. At low densities even correlated 
collisions happen at large distances, i.e. in the order of a mean free path length apart from 
each other. Therefore their distribution of collision angles and hence their contribution to the 
average ( ) c , to the leading order in density remains the same as for uncorrelated collisions. 
We should then expect that at low densities the main effect of the correlated collisions on 
the Lyapunov exponents should be due to a change of the collision frequency v as a result 
of correlated collisions taking place in the presence of the field. If the latter changes from u 
to vq + Sis, then Eq. ([15]) predicts a change of the positive Lyapunov exponent of magnitude 



To obtain this result we have used the fact that the equilibrium, low density Lyapunov 




(77) 



exponent, Eq. [15p can be written in the form 




(78) 
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where vq = 2nav. In order to understand why and how the thermostatted field changes 

the collision frequency we first recall that in equilibrium the collision frequency can be 

obtained simply by using the uniformity of the equilibrium distribution for the point particle 

in available phase space, with the result that v = -. One just has to consider the 

1 — Tina 2 

probability that the light particle during an infinitesimal time d t will hit one of the scatterers. 
On the other hand, at a time t after a given initial time, the probability for a collision may be 
considered to be a sum of three contributions: the collision frequency obtained by assuming 
that all collisions are uncorrelated and independent of each other, plus the probability for 
a recollision with a scatterer with which it has collided before, minus the reduction of the 
collision probability due to any collected knowledge of where no scatterers are present. In 
equilibrium the last two contributions have to cancel, as we demonstrate in the Appendix. 
In the presence of a field, however, this cancellation does not occur. This can easily be 
understood in a qualitative way following the argument that the cancellation in equilibrium 
occurs because the probability for return to the boundary of a scatterer is exactly the same 
as that for return to the boundary of a region where a scatterer could be, but in fact is not 
present (a virtual scatterer). In the presence of a field, the average velocity of the point 
particle before collision with a real scatterer will be in the direction of the field, and after 
the collision the average velocity will be anti-parallel to the field. The field will then tend to 
turn the particle around and have it move back in the direction of the scatterer. This effect 
enhances the probability of a recollision in comparison to that for an isotropic distribution 
around the scatterer. In a "virtual collision", in which the velocity does not change, the 
particle, on average, ends up downstream (i.e. in the direction of the applied field) from the 
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virtual scatterer and its recollision probability is decreased compared to that for an isotropic 
distribution. 

In the Appendix, a quantitative calculation is given based on the following two assump- 
tions: 

1) After the real or virtual collision the spatial distribution of the point particle becomes 
centered around a point at a distance of a diffusion length from the scatterer and 

2) for long times this distribution can be found by solving the diffusion equation. The 
resulting expression for 5u is 

— ln^. (79) 
A more formal, but equivalent, way to obtain this result is by extending the method 



described by Latz, van Beijeren and Dorfman [EH] for the low density distribution of time of 



free flights of the moving particle to include the contribution from ring events, so as to apply 
to a system in a thermostatted electric field. The main idea is to solve a kinetic equation for 
/(r, v, t, t), the distribution of particles at a phase point (r, v) at time t such that their last 
collision took place at a time r earlier, i.e., at time t — r. It is then easy to argue that the 
distribution of free flight times is simply the derivative of this ( "last collision" ) distribution 
with respect to t — r. We can then obtain a NESS average of the time of free flight and 
thereby calculate the field dependent collision frequency u(e) = vq + bv. Since we want to 
show that the origin of the non-analytic field dependence of both A + and A_ is rooted in 
the non- analytic field dependence of collision frequency 5i>, let us keep only the non- analytic 
field-dependent term as the leading term of the expansion of bv in the density of scatterers 
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and in the electric field strength and write 

f '2k u r ^ 



5v = 0£ ^{^L \ + ... j (so) 



I £ 

where the quantity j3 has to be determined from the NESS average of r, using the effect 

of the ring term on the NESS distribution function f(f,v,r) with k of the order of . 

To obtain this distribution function, we follow exactly the same procedure as outlined in 
Sections 4 and 5, but this time, with the variable r instead of p. Notice that, this time, 
even though the equations for corresponding fi and g^s are different, due to the difference 
in the dynamical equations for p and f during free flights and at collisions, the equations 
involving F% and GVs remain the same. The source of the non-analytic field-dependent term 
will surface again exactly from the 0(k~ 2 ) term in Eq. fl54"|). As far as this non-analytic 
field- dependent term is concerned, at the lowest order of density, the variables p and r are 
identical up to a multiplicative factor v. Both grow linearly with time in between collisions 
and are set back to (for p, almost) zero at each collision with a scatterer. One then recovers 



the corresponding non-analytic part of the NESS distribution function f27j, analogous to Eq. 
fSL^ft r) = { In ^1 (1 - 2navr) e" 2 — for r > , (81) 



1, non-analytic V ! / 4;7t 2 V 

from which (3 can be obtained to be 



< 82 > 



after which, one easily recovers the result of Eq. ( f7U| ) . 
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6 Discussion 

While much of this paper is quite technical, there are two main points that we would like 
to emphasize: (1) We have developed a method which allows an extension of the calculation 
of the Lyapunov exponents for a two-dimensional Lorentz gas to higher densities than is 
possible by means of the ELBE. (2) The logarithmic terms obtained here, while small, are 
indicators of similar logarithmic terms which are certain to appear when these calculations 
are extended to general two-dimensional gases, where all of the particles move. 

The first point allows one to contemplate a general kinetic theory for the calculation of 
of sums, at least, of all positive, or of all negative Lyapunov exponents. Such an approach 
was also indicated by Dorfman, Latz, and van Beijeren for the KS-entropy of a dilute 



gas in equilibrium, but the theory there has not yet been developed beyond the Boltzmann 
equation. The relevance of the second point can be seen if one realizes that the linear 
Navier-Stokes transport coefficients of a two-dimensional gas diverge because of long time 



tail effects, of the type discussed here ||10|| . In the general gas case therefore the logarithmic 
terms in the positive and negative Lyapunov exponents will not cancel as they do here, 
because the transport coefficients themselves should diverge as hie as e approaches zero. 
Thus the logarithmic terms obtained here should be seen as precursors of the more important 
logarithmic terms that will appear in the theory of two-dimensional gases. 

It is worth noting that the e 2 In e term results from a long range correlation in time 
between the moving particle and the scatterers that is present in both the pair correlation 
functions, G2, and g2, either of which is proportional to the square of the electric field 
strength and the inverse square of the wave number, at small wave numbers and fields. This 
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dependence is not present in the Lorentz gas in equilibrium, of course, but similar collision 
frequency arguments to those given here suggest that non-analytic terms may be present 
in the ring contributions to the positive Lyapunov exponent for trajectories on the fractal 
repeller for an open Lorentz gas. In this case the inverse system size, L' 1 , plays the role 
of — , the lower limit of k for the integration over k and one would expect to find terms of 
order L~ 2 In L in the ring term for this case. This point is currently under investigation. 

Finally we mention that neither the non-analytic terms found here, nor the excluded 
volume corrections included in the Enskog terms are able to account for the field dependence 
of the Lyapunov exponents as observed in the computer simulations by Dellago and Posch 
||. This is not unexpected since we have not been systematic in computing the density 
dependence of the coefficient of e 2 , nor have we considered higher order terms in e beyond 
order e 2 \n.e. All of the neglected terms are likely to be numerically more important than 
the ones we have kept. There is also no indication in the simulation data for the Lyapunov 
exponents of a clear presence of the interesting logarithmic term in the applied field. Such 
logarithmic terms are typically difficult to detect in simulation data, without a careful hunt 



for them p0|| . However, it may be easier to check, by means of computer simulation, the 
existence of the e 2 In e term in the collision frequency than in the Lyapunov exponents. In 
any case, we would like to emphasize that computer simulation studies of thermostatted 
systems provide very useful ways to check a number of phenomena predicted by the kinetic 
theory of moderately dense gases. 
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Appendix 

Derivation of the field-dependent collision frequency 

To derive the field dependence of the collision frequency we first approximate the probability 
of a recollision at time t as 

P ICC (t) = - [ dO f dr [ da\v-a\R(T,6,a)b f7 F (B \6) . (Al) 

2 Jo Jo Jv-a>0 

Here F {B> (9) describes the Boltzmann distribution for the velocity in the NESS. The function 
R(t, 9, a) describes the probability density for return to the circumference of a given scatterer 
in a time r just after colliding with this scatterer with scattering vector a and post-collisional 
velocity described by 9 (see Fig. 4). We have ignored a possible dependence of the collision 
frequency v on v, which would only play a role at higher orders in the density. 
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Fig. 4 : A recollision taking place after a real (solid line) or corresponding virtual (dashed 
line) collision, followed by a post-collisional excursion maintaining on average the direction 

of velocity over a persistence length l p . 

Similarly the reduction of the collision frequency at time t due to virtual recollisions can 
be estimated as 

P nc (t) = - V - [ del dr f da\v-a\R(r,9,a)F (B \9). (A2) 
2 Jo Jo Jv &>o 

In equilibrium F (B) (9) is independent of 9, so one sees immediately that both terms 

cancel, as they should. In the presence of a thermostatted field we need the explicit form of 

Fi B) (6) up to the first field-dependent order, given in Eq. fl50"D as 



FF{0) = i- 

271 



3e 

1 H cos 9 

8nav 



(A3) 



The function R(t, 9, a) for large enough r may be approximated by the product of lav 
(velocity times cross section) and the probability density for finding the point particle at the 
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position of the scatterer. For weak fields the latter may be approximated by the solution 
of a diffusion equation with a drift velocity ux in the +x-direction and an initial density 
localized at the position l p 9 with respect to the center of the scatterer. Here l p is the 
persistence length, that is, the average distance traveled by a point particle in an equilibrium 
system in the direction of its initial velocity and 9 is the unit vector in the direction of the 
velocity right after the initial collision at t — r. The persistence length may be expressed as 

/•OO 

lp— dt (v ■ v(t)). Multiplying this by the constant speed v one finds with the aid of the 
Jo 

2D 

Green-Kubo expression for the diffusion that L = in two dimensions. This assumption 

v 

for the long time distribution may be understood by imagining that the first few free flights 
after the initial collision of the particle move it over a distance in the order of a mean free 
path in the direction of its initial postcollisional velocity before it starts to diffuse by virtue 
of further collisions with scatterers. Thus for large r the distribution of the light particle 
will be centered around the point l p 9 with respect to the center of the scatterer, and the 
final point, on the surface of the scatterer, may be approximated to be at the center of the 
scatterer as well, because of low density. These arguments lead to the explicit form for the 
recollision probability given by 

_ [l p 6 +utx] 2 



R(t,0,°) = 2av € . (A4) 

Finally we need the explicit form u = for the drift velocity to leading order in the 

density, and the identity that 



II- 

z Jva>0 



V 

da \v ■ a\ bp cos 9 = cos6>. (A5) 

3 
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Then, after expanding 



l-^cos0 + 



(A6) 



we can now do all the calculations needed to obtain the leading non-analytic term in the 



field expansion of the collision frequency. We find that 



2tx 



8v = aw I d9 I dr 

lo Jo 4*Dt 



1-— cosfl 



7-1^ - 1 /, . % 3^ cos . . 

da\v-a\(b & -l)— . A7 

lbimav 



. Jv-a>0 

After performing the integrations, we recover Eq. (|79|) . Notice that the logarithm of e results 
from the cut-off on the r integration provided by the drift term in the exponential. 
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